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Abstract 

We study singularities of geodesics flows in two-dimensional generalized Finsler 
spaces (pseudo-Finsler spaces). Geodesics are defined as extremals of a certain aux¬ 
iliary functional whose non-isotropic extremals coincide with extremals of the action 
functional. This allows to consider isotropic lines as (unparametrized) geodesics. 


Introduction 

The paper is a study of singularities of geodesics flows in generalized Finsler spaces (pseudo- 
Finsler spaces). This is a natural development of an ongoing research on understanding the 
geometry of surfaces endowed with a signature varying pseudo-Riemannian metric; see mm 
mmmmm and the references therein. One of the purposes of the paper is to compare 
singularities of geodesics flows in pseudo-Finsler and pseudo-Riemannian metrics. On the 
other hand, the interest to pseudo-Finsler spaces is motivated by physical applications; see 

e-g. [HEJ. 

According to [23], by a pseudo-Finsler space we mean a manifold M, dim M = m, with 
coordinates (x{) endowed with a metric function f(xf, x t ) = F(xf, ±j)«, where F: TM — > R is 
positively homogeneous in (&i) of degree n and smooth on the complement of the zero section 
of TM (a more detailed definition is given in Sections 11.11 II.2j) . A well-known example is 
Berwald-Moor space (M, /), where f(xi, Xi) — {x\ ■ ■ ■ x n )™, n — m\ see, e.g., [BJ T3> [23] . 

The paper starts with the discussion of the notion of geodesics in Finsler and pseudo- 
Finsler spaces with n > 3 (Section [I]). Here we use the variational definitions of geodesics 
P |23]. In contrast to pseudo-Riemannian spaces (n — 2), where naturally parametrized 
geodesics of all types (including isotropic) can be defined as extremals of the action func¬ 
tional, in pseudo-Finsler spaces the similar definition is not correct for isotropic lines. The 
solution of this problem is either to exclude isotropic lines from consideration or to find a 
natural extension of the definition of geodesics. 
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In the present paper, we choose the second way. Based on a simple variational prop¬ 
erty, we define geodesics as extremals of a certain auxiliary functional whose non-isotropic 
extremals coincide with extremals of the action functional. In this direction, we have the 
following result: in the case m = 2 , all isotropic lines are (unparametrized) geodesics. 

In Section El we consider singularities of the geodesic flows in pseudo-Finsler spaces 
(M,/), where m = 2 and F is a polynomial in (xf) of degree n > 3. The main results 
are presented in Section 12.21 where we consider the case n — 3 in detail. It is proved that 
singularities of the geodesic flow are connected with the degeneracy of isotropic lines net. 
Namely, if the function F is generic, the manifold M contains two open domains M + and 
M_ separated by a curve M 0 so that at every point q e M + (resp., M_) there exist 3 (resp., 
1) different isotropic directions and isotropic lines are tangent at q G Mq. Singularities of 
the geodesic flow appear in the domain M_ and on the curve M 0 . 

Section El is devoted to a special case: pseudo-Finsler space (M,/), where M is a sur¬ 
face in n-dimensional Berwald-Moor space. The corresponding function F is a non-generic 
polynomial in {xf) of degree n. In this case, the domain M_ = 0, and singularities of the 
geodesic flow appear on the curve M 0 only. 

The author expresses deep gratitude to prof. Farid Tari (ICMC-USP, Sao Carlos, Brazil) 
for attention to the work, useful advices, and remarks. I am also grateful to the reviewer for 
many constructive comments and suggestions. 


1 Variational definition of geodesics 


1.1 Finsler spaces 


Consider a smooth (here and below, by smooth we mean C°° unless otherwise stated) man¬ 
ifold M, dim M = m , with coordinates (x*) and a function F(x;;Xj): TAf -)• R that is 
positively homogeneous of degree n > 2 in (xj) and smooth on the complement of the zero 
section of TM. 

_ _ 1 _ 

Define the function /(x,;Xj) = F n (x i ]x i ) ) which is positively homogeneous of degree 1 
in (hj). The pair (M, /) or, equivalently, ( M,F ) is a Finsler space (in the classic sense), if 
the following conditions hold: 

B. F(xi,Xi ) > 0 if \x\ | H-b |x m | ^ 0. 


-r2 


C. The Hessian of the function / with respect to (x*) is positive definite, that is, 


m 02 (~c 2 \ rn 

if Ei&i^°- 


1 , 3=1 


dijdxi 


( 1 . 1 ) 


i —1 


Here we use the letters B and C to preserve the notations from the book [23], to which we 
shall refer. The quadratic form (11.1(1 is called the fundamental tensor, and / (positive and 
smooth on the complement of the zero section of TM) is called the metric function on M. 

The metric function f(xp, ■) defines a Minkowski norm on each tangent space T X M. For 
a curve 7 : / —> M, it allows to define the length and the action functionals similarly to the 
Riemannian metrics: 

J (y \l) = J f U (xi]Xi)dt = I F n (xi]Xi)dt, ±i = ^, 

1 1 


2 


( 1 . 2 ) 




with v — 1 (length) and v = 2 (action), see, e.g., [13], 23]. As in the Riemannian case, the 
length functional Jd) is invariant with respect to reparametrizations of 7 , while J^ 2 ) is not. 

Parametrized geodesics can be defined as extremals of the action functional J^ 2 \ the 
corresponding parametrization is called natural or canonical (it coincides with the arc-length 
parametrization, where ds = /). 

Non-parametrized geodesics can be defined as extremals of any one of the functionals 
jd) and Jd) The difference between using Jd) an d j(b j s the following. In the first case we 
just simply forget the natural parametrization of the extremals of J^ 2 \ while in the second 
case the Euler-Lagrange system with the Lagrangian f(xf,Xi ) contains m — 1 independent 
equations only [23]. This reflects the fact that the length functional Jd)(y) is invariant 
with respect to reparametrizations of 7 . Using this degree of freedom and assuming that 
we deal with continuously differentiable geodesics with definite tangent directions at all 
points, one can put (at least, locally) the parameter t equal to one of the coordinates 07 , and 
consequently, reduce the Euler-Lagrange system for jd)(y). 

From now on, we use the following general notation. Let Q(xi\Xi) be a function on TM 
positively homogeneous of degree k in ( 07 ), then the formula 


defines a function on the projectivized tangent bundle PTM. For instance, put x = X\ and 
Hi = Xi, pi = dyi/dx for i — 2,..., m. This yields 


±(df\ 

dx \ dpi J 


d£ 

dyi 


f(x,Vi,Pi) 


f{Xi\Xi ) 

1 

Xi 


i — 2 ,..., m. 


(1.4) 


The passage to equation (11.41) is the standard projectivization II: TM —y PTM of the 
tangent bundle. Moreover, non-parametrized geodesics can be defined as extremals of 
with arbitrary v > 1 , on the basis of the following simple property (see, e.g., US]). 


Lemma 1 Let i/j: M —» M be a function such that if o F is C 2 -smooth on the complement of 
the zero section of TM and 'if'(s) 7 ^ 0 for all s 7 ^ 0. Then non-parametrized extremals of the 
functional 




1 ; c F{x,] i.) ilt, i: = ‘if . = L. 


(1.5) 


I 


coincide with non-parametrized extremals of where id is the identity map. 


Proof. The Euler-Lagrange equation of Tp(q) reads 


d 

dt 


o F(xi]Xi) ■ 



if' o F(xi,Xi) ■ 


dF 

dxi 


i — 1 ,..., m. 


( 1 . 6 ) 


In light of the condition F{xi]xf) > 0, every curve 7 admits the arc-length parametrization, 
that is, F(xi]Xi ) = c 7 ^ 0 along 7 . Using the arc-length parametrization in (II.5|) . after 
reducing the constant factor if' o F(xp,Xi ) = if'(c. ) in both sides of (II. 6 p . we get the Euler- 
Lagrange equation of the functional Jjd( 7 )- □ 
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Thus non-parametrized geodesics can be defined as extremals of 7 ^.( 7 ) with arbitrary 
function from Lemma [TJ in particular, of the functional 7 ^( 7 ), which is equal to J^( 7 ) 
with v = n from (j 1.2 jl . In classical Finsler spaces this extended definition of geodesics gives 
us nothing essentially new, but it can be useful for generalized Finsler spaces considered 
below. 


1.2 Generalized Finsler spaces 

A generalization of the notion of a Finsler space may be obtained if conditions B and C are 
dropped, such spaces are sometimes called special Finsler or pseudo-Finsler. Here we take 
the liberty to cite a passage from the classical book [23J (page 265): 

Again, it should be remarked that very frequently the metric function which is 
given by a homogeneous Lagrangian function of a dynamical system does not 
always satisfy conditions B and C. The singularities which may occur as a result 
of the relaxation of condition C are usually ignored, but it is well possible that 
an investigation of these singularities in connection with physical applications 
cannot be avoided and might furthermore prove to be fruitful. 

From now on, we will consider pseudo-Finsler spaces (M, F), where the function F is not 
supposed to satisfy the conditions B and C. 

The absence of the condition B brings to the existence of the isotropic hypersurface FA 
given by the equation F[xi,xf) = 0 in TM or, equivalently, F[x,yi,pf) = 0 in PTM. The 
Euler-Lagrange equation for the functional ,J ( F (7) with v < n is not defined on F, since 

_ v_ 

the derivatives of F n [xi\xf) are discontinuous on F. This explains the advantage of the 
functional J <yTL )[ 7 ) for definition of geodesics compared to J^ u \ 7 ) with v < n. 

The Euler-Lagrange equation for the functional F n > ( 7 ) reads 


^ d 2 F .. d 2 F . dF 

^ + Z_^ FFP-r- ^ j ~ fir ■ ’ ^ — 1, . . . , m, 


3 = 1 




3 =1 




or, equivalently, 


Hi^xf) — ( d 2 F . 

Xi = ■=;- : —, where H = det ( Q , Q . ), * = 1 ,..., m, 


H(xi, xf) 


diidij 


(1.7) 

( 1 . 8 ) 


and Hi are the determinants defined from (11.71) by Cramer’s rule. It is not hard to see 
that the functions H, are positively homogeneous of degree n in [xf) and H is positively 
homogeneous of degree n — 2 in [xf). 

Similarly to (11 .4 j) . the projectivization n \ TM PTM sends equation HI.8)1 to 


dyj dpj 
dx ’ dx 


dpj fdx\ 1 d fijA _ jjjX - ijjX _ 

dt \dt J x dt\x J x 3 

1 Hi -pjHi _ H,j - pjHi 

x 2 H ~ H 


i — 2,..., m, (1.9) 


where the functions IF II, are obtained from 77,77* by (II.3[) . Recall that, by formula (II. dp . 
the independent variable x is the coordinate X\. From Lemma [T] it follows that out of 
the isotropic hypersurface IF, integral curves of (jL9j) coincide with integral curves of (II.4)) . 
However, equation (II.9j) is defined in the whole space PTM. 
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Lemma 2 & is an invariant hypersurface of both equations (11.71) and (II.9j) . Moreover, in 
the case m = 2 all isotropic lines are non-parametrized extremals of the functional Jj 


Proof. After straightforward transformations, equation (ll.4|) gives a direction field which 
is parallel to 

f ~^Gw4) +f " p |> 1? (i - io) 

where p — (m — 1)(2 — -) and A", Y t are smooth functions on PTM. 

Since the vector held (ll.lOp is derived directly from the Euler-Lagrange equation (II .4j) . 
it is divergence-free in PTM except for the hypersurface & where the factor F is discon¬ 
tinuous, and the held is not defined. By Theorem 1 [S], & is an invariant hypersurface of 
the vector held 

which is obtained from (ll.lOj) by eliminating the common factor By Lemma [U integral 

curves of (II. lip coincide with integral curves of (11.91) . hence X = H , Yi = H, — p,jH\ , and 
& is an invariant hypersurface of (11.91) . We remark that every solution of (11.71) is obtained 
from a solution of (11.91) by choosing an appropriate parametrization Xi(f). Thus & is an 
invariant hypersurface of equation (II .7j) also. 

Finally, consider the case m — 2. Then dim PTM = 3 and dimJ^” = 2. The held of 
contact planes y 2 = P 2 dx cuts a direction held on 2 ^", which coincides with the restriction of 
the held (11.91) to &. Hence the projection of integral curves from the surface & C PTM to 
M are isotropic lines and non-parametrized extremals of the functional simultaneously. 


□ 


In accordance with what the previous reasonings, one can give the following definition. 


Definition 1 The projections of integral curves of equation (11.91) from PTM to M distin¬ 
guished from a point are non-parametrized geodesics in the pseudo-Finsler space (M, F). 

By Lemma [2J in the case m = 2 all isotropic lines are geodesics in the sense of the given 
definition. 

The natural parametrization of non-isotropic geodesics is defined by equation (11.61) with 
if(s) = s™ and coincides with the arc-length parametrization. In the case n > 2 the natural 
parametrization of isotropic geodesics is not defined, while in the case n = 2 it is defined by 
equation (II, 6 [) with if = id. 


2 Polynomial pseudo-Finsler metrics on 2-manifolds 

From now on, we consider the case when m — 2 and the function F is a homogeneous 

polynomial of degree n > 2 in (xj). Denote the coordinates on the manifold M by (x,y). 

_ _ 1 _ 

Consider pseudo-Finsler space with the metric function / = F" , where 

n n 

F o, y\ X, y) = a i( x , y) xU ~ l y\ F(x , y,p) = ^ a^x, y)p\ (2.1) 

Z =0 2=0 

* This statement holds true only for m = 2, while in the case m > 2 some isotropic lines are geodesics 
and some of them are not. An example is given for m = 3, n = 2 (pseudo-Riemannian metrics) in [201 . 
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( 2 . 2 ) 


the coefficients a, smoothly depend on (x,y). Then equation (11.8|) reads 


= Ei 

H ’ 



F ■ 

A XX 

F ■ ■ 

1 xy 

, Hi = 

G\ 

F ■ ■ 

1 xy 

, H 2 = 

F ■ ■ 

L XX 

G i 

~F. . 

1 xy 

F ■ ■ 

1 yy 

g 2 

F ■ ■ 

1 yy 

5 2 

~F. . 

1 xy 

g 2 


where G\ = F x - xF xx - yF xy and G 2 = F y — xF xy - yF yy . 


Lemma 3 The projectivization H: TM PTM sends equation (1 2.2 j) to 


where 


dy dp H 2 — pH\ 

^ dx' dx H 

A(x,y,p) = nFF pp - (n - l)F p 2 , 

P(x, y; p) = nF(F y - F xp - pF yv ) + (n 


P 

A’ 


1 )F P (F X + pF y ). 


(2.3) 

(2.4) 


Proof. Taking into account fl 1.9 [) . it remains to establish the equality IMMh — where 
A, P are defined in ( 12 .4p . Let us prove that H = (n — 1)A and H 2 — pHi — (n — 1 )P, i.e., 
H = x 2n ~ 4 (n — 1)A and H 2 —pH 1 = x 2n ~ 2 (n — l)P. Since both sides of the two last equalities 
can be treated as quadratic forms on ao,... ,a n with coefficients depending on x, y, it suffices 
to compare the coefficients of the monomials aiOj in the lift and right-hand sides. 

Put Eij — 1 if i 7 ^ j and | if i — j. Direct calculation shows that the coefficient of 

the monomial ayij, i + j = k, in the expression H = F xi F yy — f\ ± is 
where 


otij = (n - i)(n - 1 - i)j(j - 1 ) + (n - j)(n - 1 - j)i(i - 1 ) - 2 ij(n -i)(n- j ) = 

{n — 1) (n(k 2 — k — Aij) + 2 ij). (2.5) 

On the other hand, the coefficient of the monomial a t Oj, i + j = k, in the expression 
A = nFF pp — [n — 1 )F 2 is ^^e^p k ~ 2 , where 

Pij = n[i(i - 1 ) + j(j - 1 )) - 2 ij(n - 1 ) = n(k 2 -k- Aij) + 2 ij. ( 2 . 6 ) 

From (12.5p and ( 12 . 6 p . we have = [n — 1 that proves FI = i 2n_4 (n — 1)A. The proof 
of the equality H 2 — pH\ = x 2n ~ 2 {n — 1)P is similar. □ 


Remark 1 From formula (12.4)1 it follows that A and P are polynomials in p of degrees not 
greater than 2n — 4 and 2n — 1, respectively. For instance, 

A (x,y,p) = (2na n a n _ 2 - (n - 1 )a 2 n _ 1 )p 2n ~ i H-b 2 na 0 a 2 - (n - l)a 2 . 


For our further purposes, it is convenient to write equation (12.3p as the held 


( d d\ 

d ~r\ +PT^ + P 7^- 


AI-i-y- , , - 

\dx y dy) dp 


(2.7) 


The held (12.7p is defined in the whole space PTM including the isotropic surface &. The held 
of contact planes dy = pdx dehnes on & a direction held whose integral curves correspond 
to isotropic lines, while all remaining integral curves of the held ( 12 .7p (that do not belong 
entirely to the isotropic surface) correspond to non-isotropic geodesics. 

In accordance with Definition [fl non-parametrized geodesics in the pseudo-Finsler space 
(M, F) are the projections of integral curves of the held (12.7p from PTM to M distinguished 
from a point. Singularities of geodesics occur at the points of PTM where A(x, y; p) vanishes. 
To describe the locus of such points, we use the following lemma. 
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Lemma 4 Given a polynomial 


$(p) = Yl(P + 7i), 7. 6 *, W > 2, (2.8) 

consider the polynomial 

A(p) = n$(p)$"(p) — (n — l)<h'(p) 2 . (2.9) 

Then the following statements hold: 

(a) A = 0 if and only if 7 i = • • • = 7 n - 

(b) Suppose that 7 * 7 ^ 7 , for at least one pair i,j. Then p is a real root of the polynomial 
A if and only if p is a multiple root of the polynomial $. 

(c) If p is a double root o/4> and n > 3, then p is a double root of A. 

Proof. The implications 71 = ■ ■ ■ = 7 „ =>- A = 0 in (a) and $(p) = $'(p) = 0 =>■ 
A (p) =0 in (b) are trivial. The implication A = 0 =>- 71 = ■ ■ ■ = 7 n in (a) follows from 
(b). Indeed, assume that A = 0 holds and any two of the numbers 71 ,..., 7 n are not equal. 
By (b), A (p) = 0 implies <E>(p) = 0. Hence <f> = 0, which contradicts (j2.8j) . 

Statement (c) is also trivial: differentiating (12.911 twice, from <3>(p) = $'(p) = 0 we get 
A (p) = A'(p) = 0 and A "(p) = (2 — n)<f>"(p ) 2 7 ^ 0 if <1>"(p) 7 ^ 0 and n 7 ^ 2 . 

It remains to prove the implication A (p) = 0 and p 6 R $(p) = ^'(p) = 0 in the 
statement (b). Assume that 7 * 7 ^ 7 for at least one pair i , j and there exist p, eR such 
that <3>(p*) 7 ^ 0. Making the change of variables p h->• p — p*, without loss of generality we 
can assume that p* = 0. Then <f>(0) = 71 • • • 7 n 7 ^ 0, and 

n 

$ '(°) = $ "(°) = 2 ^cqct / $( 0 ), ai = —. 

1 ■ , ■ r )i 

i=l i<j 

Substituting the above formulae in ( I2.9j) . after straightforward transformations we get 

n 2 

A(0) = $ 2 (0) (2n ^ a i a i ~ ( n ~ !) (^2 “*) ) = -$ 2 (0)<Pn(ai, • • ■, «„), (2.10) 

i<j i= 1 

where 

n n 2 

• • • 1 ®n) ^ ^ ^ ' 

i= 1 i=l 

Let us prove that for any n > 2 the form <p n (ay,..., a n ) > 0 and <p n (ai ,..., a n ) = 0 if and 
only if «! = ••• = a n . Indeed, consider the vectors a = (07, ..., a n ) and f3 = (1,..., 1) in 
n-dimensional Euclidean space with the standard inner product. Then the Cauchy-Schwarz 
inequality ( a,/3 ) 2 < (a, ce) ((3, (3) gives the required assertion. 

By our assumption cq 7 ^ Oij for at least one pair i,j. Then <p n (ay,..., a n ) > 0, and 
equality (I2.10|l implies that A(0) = 0 <f>(0) = 0. Moreover, from (I2.9|l it follows A(p*) = 

<£>(p*) = 0 =>■ $'(p*) = 0, i.e., p* = 0 is a multiple root of $. The lemma is proved. □ 
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Remark 2 Obviously, the implication $ = $' = 0^A = 0 holds true for all polynomial <f>, 
not necessarily (12.81) . However, the inverse implication A = 0=^<f) = 0is not valid if $ has 
a complex root. The reason for this is easily ascertained: the inequality p n [a i, > 0 

is not valid if among the numbers cq some are complex. 

For example, consider the polynomial <f> = p 3 + p with the unique real root p = 0. Then 
the corresponding polynomial A = 2(3 p 2 — 1) has two real roots, none of those coincides 
with 0. Moreover, the polynomial $ = p 4 + 6 p 2 + 1 does not have real roots at all, while the 
corresponding polynomial A = 48(p 2 — l) 2 has two double roots p = ± 1. 

Lemma |4] gives a simple geometrical description of the singular locus of equation (12.3j) 
for the domain M' C M, where pseudo-Finsler space (M, F) has m isotropic lines passing 
trough every point of M', i.e., the polynomial F(p) has m real roots (taking into account the 
multiplicity and possibly including p = oo). For (x,y) G M' the function A (x,y,p) vanishes 
if and only if at least two of m isotropic lines are tangent at (x, y ) and p is the corresponding 
tangential direction. Remark that this statement is not valid for the complement of M', 
where the polynomial F(p) has complex roots. 

This question will be considered in detail for n — 3. 

2.1 Pseudo-Riemannian metrics 

By Remark Q] in the case n = 2 (pseudo-Riemannian metrics) A is a zero degree polynomial 
in p, that is, A does not depend on p. Moreover, it is easy to check that A = —4-Dpq, where 
means the discriminant of the quadratic polynomial F. Hence the locus of singularities 
of equation (12. 3 p coincides with the discriminant curve of the implicit differential equation 
F(x, y\p) = 0. It is not hard to see that the equation A(x, y) — 0 defines an invariant surface 
of the held (12.7j) filled with integral curves whose projections PTM —y M are points forming 
the discriminant curve. 

This property leads to a curious phenomenon: geodesics cannot pass through a point 
(x, y) of the discriminant curve in arbitrary tangential directions, but only in admissible 
directions p defined by the condition P(x,y,p) = 0. Generically, P(x,y,p) is a cubic poly¬ 
nomial in p and the number of admissible directions is 1 or 3. Singularities of the geodesic 
hows in pseudo-Riemannian metrics are studied in detail (the interested reader is referred to 
the papers jHl HJ [20] devoted to 2-dimensional pseudo-Riemannian metrics; similar results 
for 3-dimensional pseudo-Riemannian metrics were announced in [IB]). 

It should be remarked that the case n = 2 is exceptional from the viewpoint of Finsler 
and pseudo-Finsler geometry (n> 2). In the case n > 2, A generically depends on p and the 
notion of admissible directions does not appear. Geodesics pass through every point of M 
in all possible directions, but some directions at some points are singular. In other words, 
only points of the space PTM may have the property of being singular. 

In the rest of the paper, we consider the case n — 3 (cubic pseudo-Finsler metrics) in 
detail. 

2.2 Cubic pseudo-Finsler metrics 

Let D[p] and D[ a] be the discriminants of the cubic polynomial F(x,y;p ) and the quadratic 
polynomial A (x,y,p) in p, respectively. A direct calculation shows that D[ a] = —1 2D[p]- 

Introduce the following stratihcation of the manifold M. The open domains M + ,M_ 
are defined by the conditions Dyp\ > 0, Dy F ] < 0, respectively. Generically, M + ,M_ are 





separated by the discriminant curve Mq: D\p\ = 0, which consists of regular points (the 
cubic polynomial F has one prime root and one double root) and cusps (F has a triple root). 
By \ denote the set of all regular points of Mq, while Mqq = M 0 \M 0 i i. The discriminant 
of the quadratic polynomial A is strictly negative in M + , hence singularities of equation 
(12.3(1 occur only in M_ and Mq. Further we exclude from consideration the stratum Mop of 
dimension zero, and consider only M_ and Mop. 

In a neighborhood of every point of M \ M 0 0 the cubic polynomial F has at least one 
prime real root p*(x,y) smoothly depending on x, y. To simplify calculations, choose local 
coordinates such that the integral curves of the vector held — p* (x, y) (one of three 
families of isotropic lines) become x = const. This yields a^{x,y) = 0 and 

F = ap 2 + 2 bp + c, A = —2 (ap + b ) 2 + 6 (ac — b 2 ), = ia 2 (b 2 — ac ), 

M± = {±( 6 2 — ac) > 0, a ^ 0}, Mq\ = \b 2 — ac = 0, a ^ 0}, (2.11) 

P = 3 F[F y — F xp — pFyp ) + 2 F p [F x + pF y ). 

2.2.1 Singularities in the stratum M_ 

At every point in M_, the quadratic equation A = 0 has two prime real roots 

Pi, 2 = ——S = a/3 (ac-b 2 ), (2.12) 

CL 

and the domain M_ is filled with two transverse families of integral curves of the binary 
implicit differential equation A = 0 , which we shall call singular lines of the metric. 

Consider the curves S) C M_ defined by the equations P(x,y,Pi ) = 0, i = 1,2, where P 
is defined in ( 12 . lip . They can be also considered as the branches of the locus res(A, P ) = 0, 
where “res” means the resultant of two polynomials in p. In the space PTM , consider the 
corresponding curves 

Si = {(x,y;pi) : (x,y) E Si}, i = 1,2, 

which consist of singular points of the held (12.7(1 . 

By T q denote the family of geodesics outgoing from a point q = (x, y ). The simplest type 
of singularities of the geodesic how (codimension 0 ) is given in the following theorem. 

Theorem 1 Let q £ M_ and (g;p/ ^ Si, i — 1,2. Then there exists a unique geodesic 
passing through the point q with tangential direction pi: a semicubic parabola with the cusp 
at q. In particular, if q £ 1WL \ (Si U S 2 ), the family T g contains two semicubic parabolas 
with tangential directions pi, while geodesics with all remaining directions at q are smooth. 

Proof. If P(q-,pi) / 0, then, by the standard existence and uniqueness theorem, the held 
(12.7(1 has a unique integral curve 7 * passing through the point {q\Pi). From the conditions 
a / 0, ac — b 2 / 0 it follows that A and A p do not vanish simultaneously, see ( 12 . lip . 

Hence the curve 7 j has the hrst order tangency with the vertical direction (the vertical 
direction in the space PTM is called the p-direction, i.e., the kernel of the natural projection 
PTM —y M) at the point {q',pf), and the projection of the curve 7 j to M is a semicubic 
parabola with the cusp at q. □ 

Example 1 Let F be given by formula ( 12 . lip with a = 1 and 6 = 0. Then A = 2(3c — p 2 ) 
and P = 7c y p 2 + 4 c x p + 3 cc y . The stratums M_ and M 0j i are defined by the conditions 
c(x,y ) > 0 and c(x,y ) = 0 , respectively, and the curves Si are c x = ±2c y \/3c. 
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I. Put c(x,y ) = —x (Fig. [T] left). Then Si = 52 = 0 and the semiplane x > 0 (M+) is 
hlled with the net of isotropic lines y = ±|x 5 + const (dashed curves), while the semiplane 
x < 0 (M_) is hlled with the net of singular lines y = ±^(—x)§ + const (dotted curves). 
Cusps appear when geodesics (solid curves) are tangent to singular lines. Remark that 
geodesic pass from M_ to M + or vise versa through M 0j i (the y-axis) without singularity 
if they intersect the y-axis with any non-isotropic tangential direction p ^ 0. Otherwise, 
equation (12. 3 p has singularity. As we shall see in Section 12 . 2.21 at such points there exist a 
one-parameter family of geodesics outgoing in both domains M + and M_ with the common 
tangential direction p — 0 , and the prolongation of geodesics through Mo.i is not naturally 
defined. 

II. Put c(x, y ) = ay 2 — x with 0. Then Si ^ 0, i — 1, 2, but both curves Si do not 
pass through the origin. In a neighborhood of the origin that does not contain the curves 
Si, geodesics are presented in Fig. Q] (right). Here, for definiteness, we assume a > 0. All 
notations have the same meanings as before. 



Figure 1: Example [t] I and II (left and right, resp.). The stratum Mop (the y-axis on the 
left and the parabola on the right) is depicted with the bold solid line. Solid and dotted 
curves present geodesics and singular lines, respectively. Dashed curves present isotropic 

lines. 

The next type of singularities of the geodesic flow in the domain M_ (codimension 1) 
is connected with vanishing of the field (12.71) . This field belongs to a special class of vector 
fields whose singular points are not isolated, but form a manifold W of codimension 2 in the 
phase space. Yet such fields appear in many problems, see e.g. mmmmmmm - itis 
convenient to expressed the above condition in the following algebraic form: the germs of all 
components of the field at every singular point belong to the ideal / (in the ring of smooth 
germs) generated by two of them. 

The spectrum of the linear part of such a field (for brevity, we shall call it the spectrum 
of the field) contains only two non-zero eigenvalues Ai^, which play a prominent role in 
establishing the local normal form of the field. For instance, all components of the field (12.711 
belong to the ideal / = (A, P ), and the set of singular points W = S\ U S 2 . The eigenvalues 
Ai 2 and the corresponding eigenvectors are described by the following lemma. 
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Lemma 5 

1. The resonance Ai + A 2 = 0 holds at all points ( q\p ) G Si, i.e., A i)2 are real or pure 
imaginary numbers with opposite signs. 

2. The following conditions are equivalent: 

2.1. The eigenvalues Ai,A 2 at (g;p) G Si are not equal to zero. 

2.2. Si is a regular curve transversal to the contact plane pdx — dy = 0 at ( q\p ). 

2.3. Si is a regular curve and the direction Pi is transversal to Si at q. 

3. Generically, at almost all points ( q\p ) G S % the conditions 2.1-2.3 hold. 


Proof. Without loss of generality, suppose that q = 0 (the origin) belongs to Si and 
choose local coordinates centered at 0 that preserve the lines x = const and send integral 
curves of the vector held ^ = pi(x,y) to parallel lines y = const. The existence of such 
local coordinates follows from the general fact: if V\ and V 2 are smooth vector fields on the 
plane transversal at the point 0 , then in a neighborhood of 0 there exist local coordinates 
such that integral curves of V\ and V 2 coincide with the coordinate lines. 

Then the polynomials F, A, Drpi have the form ( 12 . 111 ) and the identities pi = 0 and 
b(x,y ) = 5(x,y) hold. Note that the first of them implies 3ac = 4 b 2 . From b 2 — ac < 0 it 
follows ac > 0. Taking into account 3ac = 4 b 2 } we conclude that none of the coefficients 
a, b , c vanishes at 0. Below, we present the proof for the stratum Si given by the equation 
3 c(c y — 2 b x ) + 4 bc x = 0 . The proof for the stratum S 2 is similar. 

Let A be the matrix of the linear part of the held (12.71) and Ai be the matrix of the 
Pfaffian system dA = 0, dP = 0, pdx — dy = 0 considered at arbitrary point (g;p) G S\, 
that is, for q G Si and p = 0: 


A = 


I j y 

0 0 
P, Pn 


Ai = 


p =0 


'A., 


0 


Ay Ap 


p p p 

1 x 1 y J p 

-1 


0 


p =0 


where 

A x | p=0 = 6 (ac x + a x c ) - 16 bb x , A p | p=0 = -4 ab, 

p x\ p=0 = c x (3c y - 2 b x ) + 4 bc xx + 3 c(c xy - 2 b xx ), (2.13) 

P p\ p =o = 4ac * ~ 6a * c + 2b ( 5c y ~ 2b x)- 

1. To prove the first statement, it suffices to show that tr A = 0. Taking into account 
the equality 3 c[c y — 2 b x ) + 4 bc x = 0 on Si and the identity b = 5 (that implies 3ac = 4 b 2 in 
a neighborhood of the origin), from (12.131) we have 

/ 46c \ 10c 

tr A = (A x + P p ) | p=Q = 10 [ac x + b(c y - 2 b x )] = 10 [ac x - ~ 4&2 ) = °- 

Hence the characteristic equation of the matrix A reads A(A 2 H- |A x |) = 0, this yields the 
equation A 2 + |Ai| = 0 for Ai i2 . 

2. Differentiating the identity 46 2 = 3ac by x, we get 8 bb x = 3 (a x c + ac x ). Using these 
identities and (I2.13p . we have 


A *| p =o = 6 ( aCa; + aa;C ) _ 1Qbbx = Q ( aCx + a * c ) _ 6 ( ac ^ + a * c ) = 

I Ai| ApD x |^_Q 1 Gab c xx Sabb x c x T \‘2ab^c x c y T cc xy 26 xx c] 

12a 2 cc xx — 3 ac x (a x c + ac x ) + 12 ab[c x c y + cc xy — 2 b xx c\. 
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Thus the condition Ap 2 7 ^ 0 is equivalent to |Ai| 7 ^ 0 that, in turn, is equivalent to the 
condition 2 . 2 . 

On the other hand, the curve Si is tangent to the direction p — 0 at the point q — 0 
if and only if [3 c(c y — 2 b x ) + 4 bc x ]' x = 0. Taking into account the equalities 4 b 2 = 3 ac and 
8 bb x = 3 (a x c + ac x ), we get 


r / , , , Ab 2 c TT — 2 bb r c r 

[3 c(c y - 2 b x ) + 4 bc x \ x = 3 (c x c y + cc xy ) - 66 ra c H---- 


3 ( 0 ^ 0 ^ T cc X y^j Qb xx c + 


12acc X3; — 3(a x c + ac x )c x |Ai| 


46 


4a6 


This proves that Ap 2 7 ^ 0 is equivalent to the condition 2.3. 

3. Generically, at almost all points (q;p) G Si the determinants 


A, Ay 

P x Py ’ 


I Ai| = 


A x A % 
P, P n 


are not equal to zero. Hence S\ and S\ are regular curves and, moreover, the conditions 
2.1-2.3 hold. n 

Theorem 2 Let {q]pf) G Si be a generic singular point of the field 02.71) . Then the germ 
02 .7|) at ( q\pi ) is smoothly orbitally equivalent to 

^ ~ ^ + tf Al ’ 2GM \°’ ( 2 ‘ 14 ) 

~ + + r] 2 " > ~d( : ^ Al ’ 2Gl \ 0 ’ ( 2 - 15 ) 

where M, I are real and imaginary axes, respectively. 

In the first case there exist two geodesics passing through the point q G Si with the 
tangential direction pi, both of them smooth. In the second case, there are no geodesics 
passing through the point q G Si with the tangential direction p^. 


Proof. Since {q\pf) £ Si is a generic singular point, |Ai| 7 ^ 0. By Lemma El the 
eigenvalues Ai ;2 0 and S t is a regular curve consisting of singular points of the field 02.71) . 
The linear part of the germ 02.7|) at (q]Pi) (and every singular point sufficiently close to 
( q\pi )) is orbitally equivalent to the linear part of the field fj2.14j) or 02.15p if |A x | < 0 or 
|Ai| > 0, respectively. Here we use the following terminology: two vector fields are called 
orbitally smoothly (resp. topologically) equivalent, if there exists a diffeomorphism (resp. 
homeomorphism) that conjugates their integral curves, i.e., orbits of their phase flowsJU 
Recall that the field 02.71) belongs to the class of vector fields whose singular points are not 
isolated, but form a manifold W of codimension 2 in the phase space (in our case, W = Sf). 
Local normal forms of such fields were studied by many authors. I 11 |21j (Appendix A), we 
present a brief survey of such results, which covers all cases with Re Ai i2 7 ^ 0. This condition 
is equivalent to the assumption that W — S t is the local center manifold, and consequently, 

t It slightly differs from the generally accepted definition of the orbital equivalence, where coincidence 
of the orientation of integral curves is also required. Our definition is naturally related to directions fields, 
whose integral curves do not have an orientation a priori. 
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the phase portrait of (12.711 has a simple topological structure (we shall discuss it later on, in 
the proof of Theorem [I]). For instance, in the case A 12 G R \ 0, the germ (12.71) with generic 
quadratic part is smoothly orbitally equivalent to (12.14j) . This result belongs to Roussarie 
1221 . The genericity condition is determined explicitly in [ 2T] (Theorem 5.7). The case 
A 12 € I \ 0 is more complicated. However, in [3] (Chapter 2, Section 1.2) and [12] it claims 
that in this case the germ (12.71) with generic quadratic part is smoothly orbitally equivalent 
to (12H5D . 

Remark that the diffeomorphism that brings the germ (12.7|) to the normal form (I2.14p or 
(12.15)) does not give a normal form of equation (12.3[) . since it does not preserve the contact 
structure dy = pdx. However, we need not a normal form of (12.3[) . 

To prove the last statement of the theorem, we only need to consider the possible mutual 
relationships between the phase portrait of the germ (| 2 .7[) and the ( x , p)-plane in the space 
PTM. Geodesics are obtained from those integral curves of the held (12.7j) whose projection 
on the ( x , p)-plane are distinguished from points. Moreover, isotropic geodesics correspond 
to those integral curves that belong to the isotropic hypersurface & (by Lemma [2j & is an 
invariant hypersurface of the held (12.7)) ). 

We consider the real and imaginary cases separately. 

The real case. The held ()2. 1 4j) has the hrst integral £ 77 . The invariant foliation £77 = const 
contains only two leaves £ = 0 and 77 = 0 that pass through singular points of the held, while 
all remaining invariant leaves are hyperbolic cylinders £77 = const 7 ^ 0 , which do not intersect 
the set of singular points. It is easy to see that for every singular point of the held (12.14)) 
there are only two integral curves passing through this point: the straight lines parallel to 
the £-axis and the 77 -axis, respectively. 

We prove now that the eigenvectors with the eigenvalues Ai )2 7 ^ 0 are not vertical. Let 
e be an eigenvector of the matrix A with A*. Then Ae = A^e, and e = ad x + (3d p , where 
(A„| p —At)di + Ap| ^_ 0 /3 = 0 . If the eigenvector e is vertical, i.e., a = 0 , /3 7 ^ 0 , this equality 
yields A p | 0 = 0. From ()2.13j) . we have a(0) = 0 or 6(0) = 0. This contradicts to the fact 
(established in the proof of Lemma [5]) that none of the coefficients a, h, c vanishes at 0. 

From the considerations above, it follows that the held (12. 7 j) has only two integral curves 
passing through the given point (g; Pi ), both of them smooth and have non-vertical tangential 
directions. Projecting these integral curves from PTM to M, we get two smooth geodesics 
passing through the point q G S) with the tangential direction p*; see Fig. [2] (left). 

The imaginary case. The held (12.151) has the hrst integral £ 2 + 77 2 . The invariant foliation 
£ 2 + 77 s = const contains a one-dimensional degenerate leaf £ = 77 = 0 , which consists 
of singular points of the held (12.15)1 and one-parameter family of two-dimensional leaves 
(elliptic cylinders £ 2 + rf = const 7 ^ 0 ), which do not intersect the set of singular points. 
The elliptic cylinders are filled with helix-like integral curves, whose projections to M have 
cusps; see Fig. [2] (right). 

To complete the proof, observe that in both real and imaginary cases the curve Si itself 
is not a geodesic, since S t is transversal to the contact plane pdx — dy = 0 (statement 2.2 in 
Lemma [3]). Consequently, Si is not a lift of a curve on M. □ 

Example 2 Let F be given by formula (12.111) with a = 1, b = 0, c(x,y) = ay 2 — x. Then 
A = 2(3c — p 2 ), P = Vlayp 2 — 4p + 6 apc, and the curves Si are the connected components of 
the graph x = ay 2 — -^/ a 2 y 2 lying the the upper and lower semiplanes. A straightforward 
calculation shows that p*| s = 12 ay (ay 2 — x), hence the direction p l is tangent to the curve 
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Figure 2: The phase portraits of the field (12.71) with the normal form (I2.14|) or (12.15|) (left 
and right, resp.) and the projections of its integral curves to M. Dashed curves present S t 
(up) and Si (down). 

Si at x = %/y/o. only. By the statement 2.3 in Lemma 0 the eigenvalues A 12 ^ 0 at all 
points of Si if a < 0 and at all points of S t with x ^ %/y/oi if ol > 0. 

In Fig. [3] (left and center) we present geodesics in the case a > 0. Here both real and 
imaginary eigenvalues exist. The parts of Si with real (imaginary) eigenvalues Ai ; 2 7 ^ 0 are 
presented as short-dashed (resp., long-dashed) lines. The dots represent the points of Si 
with x = | ^/y/a, where A i;2 = 0. In Fig. [3] (right) we present geodesics in the case a < 0. 
Here only real eigenvalues exist, and the curves Si are presented as short-dashed lines. The 
dots represent the points where geodesics intersect the curves Si with the singular tangential 
direction p t . 



Figure 3: Example [21 F = p 2 + c, where c = ay 2 — x with a > 0 (left, center) and 
a < 0 (right). The solid lines present geodesics. M 01 is depicted with bold solid line, S t are 
depicted with short-dashed (long-dashed) lines if Ai )2 £ M \ 0 (Ai >2 6 I\0, respectively). 


14 


















2.2.2 Singularities in the stratum M 0 p 

In this section as before, we proceed in the local coordinates where F, A and D\m have the 
form (12.11(1 . At every point q G M 0> 1 the polynomial F has the double root po = —b/a. It is 
easy to see that po is also a double root of the polynomial A at q (it follows from Lemma [4] 
also). Thus (q;po), q G Mop, are singular points of both implicit differential equations F = 0 
and A = 0. From (12.4p it follows that P(q-,p 0 ) = 0, hence ( q;p 0 ), q G M 0 p, are singular 
points of the field (12.711 . 

Further we restrict ourselves to generic points q G M 0 p where M 0 p is a regular curve 
and the isotropic direction po is transverse to M 0 p. Then both implicit differential equations 
F — 0 and A = 0 have Cibrario normal form at such point and their integral curves are 
semicubic parabolas lying on opposite sides of M 01 (in the domains M + and M_, resp.) as 
it is presented in Fig. [Q 


Theorem 3 Suppose that the isotropic direction po is transverse to M 0) i at q G Mop. Then 
the germ (12.7(1 at its singular point (g; po) is smoothly orbitally equivalent to 

d d d . . 

3i r ( +2 % +0 rc (216) 

and to p 0 corresponds a one-parameter family of geodesics outgoing from q into M + and M_. 
There exist smooth local coordinates centered at q such that this family is 

x = a\r)\% + p 2 + aX Q (rj), y = ar)\r)\% + erf + aY a (p), e ^ 0, (2.17) 

where X a (p ) = o(|r/|f) and Y a {rj) = o(|//|i) are C 2 -smooth functions. 

Here a > 0 (a < 0) corresponds to non-isotropic geodesics outgoing from q in M + (resp., 
M_), while a = 0 gives the isotropic geodesic, a semicubic parabola lying in M + . The limit 
case a —)■ oo corresponds to a unique smooth geodesic passing through q with the direction 
Po. In a neighborhood of q, every non-isotropic geodesics outgoing from q in M + belongs to 
the curvilinear tongue-like sector bounded by the branches of the isotropic geodesic as it is 
presented in Figure\f\ (left). 


Proof. Without loss of generality, suppose that q = 0 (the origin of the (x, y)-plane) and 
choose local coordinates centered at 0 that preserve the lines x = const and give b(x,y ) = 0 . 
It can be done using appropriate change of variables y hG yu(x,y), w(0) ^ 0. Remark that 
(unlike Lemma [3]) we do not have the identity p\ = 0 nor p 2 = 0 in a neighborhood of 0. 
Moreover, it is impossible to get any these identities using smooth change of variables, since 
the integral curves of the implicit equation A = 0 with roots Pi,P 2 have cups on Mop. 

In the local coordinates chosen above, we have 

F = ap 2 + c, A = -2 (ap) 2 + 6 ac, D [F] = -4 a 3 c, ^ 

P = aa y p A — 2aa x p 3 + (7 ac y — 3a y c)p 2 + (4 ac x — 6 a x c)p + 3 cc y . 

The curve M 0 p is given by c(x,y ) = 0 and the direction p 0 = 0 at every q G M 0 p. Hence 
the condition “the direction p 0 is transverse to M 01 at 0” is equivalent to c^O) ^ 0. 

Substituting A and P from (|2. 18jl into (12.7(1 . one can find that the spectrum of the 
field (12.7(1 at every point (g;0), q G M 0 p, is (Ai,A 2 , 0 ), where Ai = 6 ac x , \ 2 = 4 ac x . A 
straightforward computation shows that the corresponding eigenvectors are 


c i 2ad x T 3c y 0p , e2 Co c. y d x e x d y 


(2.19) 
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Figure 4: Illustrations to Theorem [3] (left) and Theorem [4] (right). The stratum Mo,i is 
depicted with the bold solid line. Solid and dashed curves present non-isotropic and isotropic 
geodesics, respectively. 

Note that Ai : A 2 = 3 : 2 at all points (g; 0), q G M 0jl , the pair (Ai,A 2 ) is non-resonant 
and belongs to the Poincare domain. Therefore, the germ (12.71) at 0 is smoothly orbitally 
equivalent to the linear held (12.161) (Theorem 5.5 in [2Tj). Moreover, comparing (12.71) and 
(I2.16|) . one can see that the conjugating diffeomorphism ( x,y,p ) i —> can be chosen 

in the form 

x = 2a£ + c y ( + P = 3c„f + p + / 2 (f, p, C), y = -c x ( + / 3 (f, p, C), ( 2 . 20 ) 

where a, c x , c y are evaluated at 0 and /) G DJI 1 (DJl k , k > 0, is the ideal of k-fi .at functions in 
the ring of smooth functions). 

The held (j2.16[) has the invariant foliation ( = const. The set of integral curves of (12.16[) 
passing through the origin consists of the £-axis and the one-parameter family 

{£ — a \v\K C — 0}, «GR, ( 2 . 21 ) 

tending to the ^-axis as a —> 00 . Consider the possible mutual relationships between the 
phase portrait of the germ (12.71) at 0 and the (x, ?/)-plane in the space PTM using the 
eigenvectors (j2.19j) . The integral curve of the held (12.7|) corresponding to the £-axis in (12.161) 
has a non-vertical tangential direction at 0 (the eigenvector ei), hence its projection to the 
(x, ?/)-plane is a smooth geodesic. On the contrary, the family (I2.21j) gives a family of integral 
curves of (12.71) with vertical tangential direction at 0 (the eigenvector e 2 ). The projections 
of these curves to the (x,y)-plane have singularity at 0 . 

To establish the character of the singularity, substitute (I2.2ip in (I2.20|) . This yields 

3 _ _ 

x = 2aa\ri\i + f 1)a (rj), p = rj + f 2 , a (v), 

where fi, a (v) — fi( a \ r )\^:' r h 0)- Observe that the functions /i )Q = o(|r/|^) and / 2)Q . = o(r]) are 
C 2 and C 1 , resp. Denote the sign of q by s(rj), then we have the equation 

dy = pdx = (p + / 2 , Q (??)) (3aa\r]\h(r]) + f[ i0t (rj)) dp = (3aa\r]\% + g a {p)) dp, 

where g a = o(|r/|i) is C' 1 -smooth. Integrating, we get y = ^aa\p\is(p) +h a (p) = ^aap\p\% + 
h a (p), where h a = o(M^) is C 2 -smooth. The scaling y 1 —y | y, a 1 —> ±2 aa yields 

x = a\p\* + X a (p), y = ap\p\i +Y a (p), (2.22) 
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where X a = o(|r/|i) and Y a = o(|r/|i) are C' 2 -smooth, and a > 0 (a < 0) corresponds to the 
domain M + (M_, resp.). The asymptotic formula (I2.22p makes sense for all real a^O. 

In order to take care of the omitted case a = 0, recall that the isotropic surface & is an 
invariant surface of the field (I2.7|) (Lemma |2J) and contains its singular points (Lemma 0J) . 
Hence in the normal coordinates (£,??, () the surface & contains the ((-axis and intersects 
every invariant leaf ( = const by a certain integral curve of (12.161) . For instance, & intersects 
the leaf ( = 0 by an integral curve of the family (12.2ip . which corresponds to an isotropic 
geodesic passing through 0. 

On the other hand, we know that the implicit differential equation F = ap 2 + c = 0, 
which described the isotropic lines in ( M,F ), has Cibrario normal form at 0. Hence there 
exists a unique isotropic geodesic passing through 0, the semicubic parabola 

x = V 2 , V = V* N (v), N( 0) ^ 0, (2.23) 


lying in the domain M + . 

From the uniqueness of the isotropic geodesic passing through 0 it follows that the lift of 
(I2.23|) is the curve of the family (12.22) ) with a = 0 and Yq(t)) = r] 3 N(r]). Using the represen¬ 
tation N(r /) = Ni(rj 2 ) + r}N 2 (rj 2 ) and the change of variables y t->- Ni(0)(y — x 2 N 2 (x)) / Ni(x), 
we get N(rj) = N( 0). It is not hard to check that the number £ = iV(0) is a mutual invariant 
of the curves (12.22j) and fl2.23[) . q 

In Example [H we considered F = p 2 + c with c = — x and c = ay 2 — x. In both cases the 
isotropic direction po = 0 is transverse to the curve Mo.i given by the equation x = 0 and 
x = ay 2 , respectively, and the conditions of Theorem [3] hold true. 

Example 3 We consider here the case F = p 2 — x in more detail. The held (12.71) is 

- 2(3 x+p2) {ik + %)- 4p w (2 - 24) 

It is easy to check that the isotropic surface & given by p 2 = x is an invariant surface of 
the held (12.241) and the unique isotropic line passing through 0 is given by x = p 2 , y = |p 3 . 
Integrating the equation dp/dx = 2p/(3x + p 2 ), we get the family x = a\p\% + p 2 , where 
a is the constant of integration, and a single integral curve p — 0, which gives the smooth 
non-isotropic geodesic y = 0. 

Integrating the relation dy = pdx = p(\ot\p\* s(p) + 2 p)dp = (|a|p|a + 2 p 2 )dp, we get 
y = \ap\p\ ^ + |p 3 + ci, where c\ is the second constant of integration. The family of geodesics 
outgoing from q — 0 is characterized by C\ = 0. The scaling y —>• |?/ brings this family to 
the form (12.171) with X a = Y a = 0 and rj = p: 

x = a\r)\ 2 + rf, y = arj\rj\^ + ^r/ 3 . (2.25) 

For a = 0, formula (12.25)) gives the isotropic geodesic. For a —> oo, the curves (12.25)) tend 
to the smooth geodesic y = 0. 

Remark 3 Theorem [3] shows that the extension of geodesics through the curve Mo.i is not 
uniquely defined. Indeed, all geodesics of the family (12.17)) have the same tangential direction 
at q G Mo,i and almost all of then have singularity of the same type at q. So, a curve given 
by formula (12.171) with any a ^ 0 does not have any advantages in comparison with the 
curve consisting of two bows (12.171) with aq if (x,y) G M + and a 2 if (x,y) G M_. 
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3 Surfaces in Berwald—Moor spaces 

Consider the space M™, n > 3, with the coordinates (aq,..., x n ) equipped with Berwald-Moor 

metric ds = (dx i • • • dx n )™, and a smooth two-dimensional surface M C parametrized 

by Xi = fi(x,y), i = The Berwald-Moor metric of the ambient space defines 

_ _ _ 1 _ 

two-dimensional pseudo-Finsler space ( M,F ) with the metric function / = F n , where 

F(x,y;x,y) = f[{fix(x,y)x + f iy (x,y)y), = f iy = (3.1) 

i= i y 

and n families of isotropic lines 

fi(x,y) = const, i = l,...,n. (3.2) 

Given q £ M, the isotropic direction p is called simple ( double or multiple ) if there exist 
only one ( only two or more than one, resp.) isotropic lines (13.2(1 passing through q with given 
direction p. By Lemma SI singularities of the geodesic flow occur at the points q £ M that 
have at least one multiple isotropic direction. 

Remark 4 In the case n — 3 we have a cubic pseudo-Finsler space ( M,F ). But unlike 
Section [2721 the function F given by (13. 1 p is not generic. The corresponding cubic polynomial 
F(x,y,p ) at every point q £ M has n real roots (taking into account the multiplicity and 
including the root p = oo), and M_ = 0. Hence singularities of geodesic appear only at the 
points where at least two of three isotropic lines (13. 2 p are tangent. Here the stratum M 0j i 
consists of the points where two isotropic lines are tangent (the double isotropic direction) 
and the third one is transversal to them (the simple isotropic direction). 

From now on, we assume that the functions fi(x,y) have non-degenerate differentials 
and every point q £ M may have simple or double isotropic directions only (the number of 
double isotropic directions can vary from 0 to [^]). Moreover, assume that the tangency of 
isotropic lines with double isotropic directions has first order. Consider geodesics passing 
through a point q with a double isotropic direction p 0 satisfying the above conditions. 

Without loss of generality assume that q — 0 (the origin in the (x, ?/)-plane) and p 0 
corresponds to the isotropic lines fi(x,y) = 0 and f 2 (x,y) = 0, where f 2y ( 0) ^ 0 . Making 
the change of variable y K > f 2 (x, y), we transform the metric function (13.1|) into a similar one 
with f 2 (x,y ) = y and fi x (0) = 0, fi y (0) ^ 0, fi xx { 0) ^ 0. The double isotropic direction p 0 
becomes p — 0 and, moreover, in a neighborhood of 0 , p = 0 is the double isotropic direction 
at all points q £ M 01 = {fi x (x, y) = 0 }. 

By the condition f\ xx ( 0) ^ 0, M 0j i is a smooth curve transversal to the x-axis. Making 
the change of variable x i —> f\ x {x,y), we transform M 0i i into x = 0 and the metric function 
(13.1(1 into a similar one with j\ x = xa(x, y), f\ y = b(x, y), y) = y, where a, b are smooth 
functions non-vanishing at 0. So, we get 

n 

F(x,y;p) = p(ax+ bp)G, G(x,y;p) = IJ(/i® + fi y p), G(0,0;0)^0. (3.3) 

i=3 

Substituting (13.31) in (12.41) . we get 

A = ((1 — n)(ax) 2 + 2(2 — n)abxp + 2(2 — n)(bp) 2 ) G 2 + A 0 , 

P = ap[(n — 2 )bp — ax + P Q ) G 2 , 
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where A 0 G (x 3 ,x 2 p,xp 2 } p 3 ) and Po G (x 2 ,xp,p 2 ) (both ideals are in the ring of smooth 
functions on x,y,p). Formula (I3.4j) shows that all components of the field (12.71) vanish on 
the line {x = p = 0} and the spectrum of (I2.7|) at any point of this line has three zero 
eigenvalues. This does not allow to establish a normal form similarly to Theorem |3j 
To overcome this problem, consider the blowing up 

B\ (x, y, u) (->■ (x, y,p), p = xu, n G IP 1 = M U oo. (3.5) 

The mapping B is one-to-one except on the plane fl = {(x, y, m) : x = 0}, whose image is the 
line £>(hl) = {(x,y,y ): x = p = 0}. The mapping B is a local diffeomorphism at all points 
of the (x,y,u)-space except fl. It has an inverse defined on R 3 \ £>(II) given by 

B^ l {x,y,p) = (x,y,~). 

\ CC / 

Observe that there is no geodesic that coincides with the line B(J\). A straightforward 
calculation shows that the held (12.71) corresponds to a smooth held in the (x, y, u)-space 
(away of fl) that, after dividing by the common factor xG 2 , is 

X ( A + • • •) (J^ + xu §^j + < n ~ 2 ) ( 2 (M 2 + 3abu + ° 2 + • • •) ^ ( 3 - 6 ) 

where 

A(x, y, u) — (1 — n)a 2 + 2(2 — n)(abu + ( bu ) 2 ), 
here and below the dots mean terms that belong to the ideal (x). 

Remark 5 There exists e > 0 such that A(x,y,u ) ^ 0 for all u if |x| + \y\ < e. Indeed, 
consider A(x, y , u) as a quadratic polynomial on u with the discriminant 

D = (2 — n) 2 (ab) 2 — 2(1 — n)(2 — n)(ab) 2 = n( 2 — n)(ab) 2 , 

which is strictly negative if x, y are sufficiently close to zero. 


Dividing the held (I3.6|) by (A +...), we get 


( d d\ 

x — + xu— + u[n 
\ox oy J 


2 )(U+.. 



U(x,y,u) 


2 (bu) 2 + 3 abu + a 2 
A(x,y,u) 


(3.7) 


Remark 6 The plane x = 0 is invariant for the helds (13.6|) and (13.70 . Moreover, it is filled 
with vertical (i.e., parallel to the ^-direction) straight integral lines of these helds, whose 
projections to the (x, y)-plane along the u-axis are points on the y-axis. 


Lemma 6 Geodesics can pass through a point q G M lying on the y-axis with the direction 
p = 0 only with the following admissible values u : 

«o = 0, ui = —a/26, U 2 = —a/ 6 . (3.8) 

Proof. By the standard existence and uniqueness theorem, for every point (x, y, u) such 
that x = 0 and U (x, y, u) ^ 0 there exists a unique integral curve of the held (13.71) passing 
through this point. By Remark | 6 l it is a vertical straight line, whose projection to the 
(x, y)-plane is a point on the y-axis. Hence geodesics can pass through a point q G M lying 
on the y-axis with the direction p = 0 only with u = 0 or u such that 2 (6m) 2 + 3abu + a 2 = 0. 
This gives the three values in (13.81) . □ 
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Lemma 7 The set of singular points of the field, (IX7D consists of three mutually disjoint 
curves 

Wf = {(x,y,u):x = Q,u = u^y)}, i = 0,1, 2 . 

On every curve Wf, the linear part of the field (13.71) has the constant spectrum (1,A,0), 
where A = rj —^ if i = 1 or A = if i G {0, 2}. In both cases, d u is the eigenvector with A. 

Proof. The first statement is trivial. All other statements are by direct calculations. □ 

Theorem 4 Suppose that the functions fi(x,y ) have non-degenerate differentials and p 0 is 
a double isotropic direction at q G M such that the corresponding isotropic lines have first 
order of tangency at q. Then the field (I3TD at its singular points (q; uf) has local orbital 
normal forms indicated in Tableland to p 0 corresponds a one-parameter family of geodesics 
outgoing from q. There exist smooth local coordinates centered at q such that this family 
consists of C 2 -smooth non-isotropic geodesics 

y = x 2 + Y(x, o|x| A ), Y(x,a\x\ x ) = o(x 2 ), \ = a G R, (3.9) 

where Yf,-) is a C 00 -smooth function, together with two C°° -smooth isotropic geodesics 

y = 0 and y = 2x 2 + o(x 2 ). (3.10) 

In a neighborhood of q, every geodesic of the family (13. 9 1) belongs to the curvilinear tongue-like 
sector bounded by the curves (I3.10|) as it is presented in Figure [7] (right). 



Orbital normal form 


topological 

C' 00 -smooth 

Wf 


^4 + A^ + 0 |, where A = ^ 

Wf 

and 

Wf 

£JL - T /— + 0 -^ 

1 drj ^ U df 

£(» -1 + »i(p, C»i - n(n - 2 + <))& + p'Hfi, Owe 

where p = 

(n - l)f| - (n - 2 )i?|; + if *( 0 , 0 ) # 0 . 


Table 1: Local orbital normal forms of the field (13.71) . 


Proof. Choose local coordinates so that q G M is the origin and consider the field (13.7|) 
in a neighborhood of its singular points ( 0 , 0 , Mj), i — 0 , 1 , 2 , where Ui are given by formula 
(13.8[) . By Lemma 0 in all singular points the condition ReAi^ 0 holds, and every curve 
Wf, i = 0,1, 2, is the center manifolds of this field. Moreover, there exist also 2-dimensional 
unstable manifold if % = 1 and the pair of 1 -dimensional stable and unstable manifolds if 
i — 0,2. Hence all topological normal forms in Table [0 trivially follow from the reduction 
principle [21 SI E] - 

Indeed, the reduction principle asserts that the germ (13.71) is orbitally topologically equiv¬ 
alent to the direct product of the standard 2 -dimensional node (if i = 1 ) or saddle (if i = 0 , 2 ) 
and the restriction of the held to the center manifold Wf. Since the restriction of the held 
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(j3.7|) to every center manifold Wf, i = 0,1, 2 , is identically zero, this gives us the topological 
normal forms in Table [I] 

Establish now the smooth normal forms in the cases i — 1 and i = 0,2 separately. 

The case i = 1. By Lemma [TJ the linear part of the held (13.711 at any point on has 
spectrum (l,A,0) with A = —Then Theorem 5.5 in m asserts that the germ (13.711 at 
any point on Wf is orbitally C'°°-smoothly equivalent to 

(£ + T(C)V 1/X )-^ + + 0|, (3.11) 

where </?(£) = 0 if the number 1/A is not integer (non-resonant case). 

Assume 1/A is integer and prove that <p(() = 0 iff for every point cu* G the held (I3.7|l 
has a C^-smooth integral curve passing through tu* with the vertical tangential direction d u . 
By Remark | 6 j such integral curve exists (the vertical straight lines), hence this establishes 
the equality </?( Q = 0 in the remaining cases 1/A = 3 (n = 3) and 1 /A = 2 {n — 4). 

For this, note that the held (13. lip has the invariant foliation = const. Every invariant 
leaf contains a single integral curve r] = 0 corresponding to eigendirection with the eigenvalue 
1 and one-parameter family of integral curves 

f = 77 1/A (a + ^(C)ln|77|), a G M, (3.12) 

corresponding to eigendirection with the eigenvalue A. All curves (13. 12j) are C 1 1 -smooth 
(but not C 1//A -smooth at zero) if ^ 0 and C°°-smooth if <p(() = 0. Without loss of 
generality, assume that the point w* G Wf in the (x,y,u)- space corresponds to (0,0,</*) in 
the (£,??, C)-space. The equality (^(C*) = 0 is equivalent to the existence of at least one C°°- 
smooth integral curve of the held (13. lip with tangential direction d v lying on the invariant 
leaf {C = £*}• To complete the proof, remark that the eigendirection d v of (13. lip corresponds 
to the eigendirection d u of (13.7p . 

The cases i = 0,2. By Lemma [7] the linear part of the held (13.711 at all points on the 
curves Wq and hL 2 c has spectrum (Ai,A 2 , 0 ), where Ai = 1 and A 2 = This gives the 

resonance 

/iAi T z/A 2 = 0 (3.13) 

with the resonant monomial p = where we set y = n — 2 and v — n — 1. Everything 
that we say below is true as well for arbitrary relatively prime /i, v G N. 

The resonance (13. 13p does not allow to get a normal form with one identically zero 
component (as we have in the case i = 1) even in the f ini te-smooth category, see the discussion 
in [8j (Section 3.2). Moreover, (13. 13p generates two infinite series of resonances 

(1 + Z/x) Ai + lv A 2 = Ai, //iAi + (1 + l is) A 2 = A 2 , l = 1 , 2 ,..., 


and consequently, infinite number of resonant monomials in the corresponding (orbital) 
Poincare-Dulae normal form: 


£(z/ + $i(p,C)) 


d_ 


r/(/i + $ 2 (p,C)) 


d_ 

dp 


+ P^(pX) 


d_ 

d? 


(3.14) 


see e.g. [H] (Section 3.2) or [T7] (Section 5). 

Moreover, if in addition, $ / 0 at a point cu*, the germ (13. 14p at a;* is smoothly orbitally 
equivalent to 


J d d 


(3.15) 
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The normal form (13. 15 ft was firstly established by Roussarie [22] in the partial case p = v = 1 
in C'°°-smooth category. For arbitrary integers p, u, the proof (in finite-smooth category) can 
be found in [T7] (Section 5). Combining the methods from [52] and [17], one can establish 
the normal form (13. 1 51) with arbitrary /i, v in C^-smooth category also. 

Completion of the proof. Integral curves of the field (13.71) passing through (0,0, Mi) 
correspond to integral curves of the field lying on the invariant leaf {£ = 0 }: a 

single curve that coincides with the r ]-axis and one-parameter family {77 = a|£| A , ( = 0 }, 
a G K. Comparing the germ (13.71) at (0, 0, uf) with its normal form fjh _|_ one can see 

that the conjugating diffeomorphism (x,y,u) 1 —» (£, r/, £) can be chosen in the form 

x = Z, u = u 1 + ri + c 1 z + c2( + <p{£,'n,C), y = t + W^vX), (3.16) 

where Wl k , k > 0, is the ideal of k -flat functions in the ring of smooth functions. Substituting 
£ = C — 0 in (13.16ft and taking into account p = ux, we get x = p = 0 . Hence the //-axis 
does not correspond to a geodesic. 

Substituting p = ct|£| A and ( = 0 in (I3.16|) and taking into account x = £ and p = ux, 
we get p = x(u\ + f(x,a\x\ x )) with a certain smooth function / e 9H°. This gives the 
relation dy = pdx = x(u± + f(x, a\x\ x ))dx, where xf(x,a\x\ x ) = o(x) is a C 1 -smooth 
function. Integrating, we get y = ^fx 2 + Y(x,a\x\ x ). Here Y(-, •) is a smooth function and 
Y{x,a\x\ x ) = o(x 2 ) is C 2 -smooth. After the scaling y —> 2y/u\, we get the family (13.91) . 

The topological and smooth orbital normal forms in Table |T] show that the field (13.71) 
has only two integral curves passing through its singular point ( 0 , 0 ,Wj), where i = 0 or 
2. Moreover, one of these integral curves is straight vertical line, whose projection to the 
(x, 7 /)-plane is a point (see Remark [ 6 ] and Lemma [7]). Another integral curve has non-vertical 
tangential direction at ( 0 , 0 , uf), hence its projection to the (x, r/)-plane is regular. 

Thus every of the admissible values u 0 and u 2 gives a smooth geodesic passing through the 
point q with tangential direction p — 0. It is not hard to see that these geodesics are isotropic 
lines, which are solutions of differential equations p = 0 and ax + bp = 0 , respectively (see 
formula (13.31) 1. Taking into account (13.8ft . after the scaling y —> 2y/u\ we get (13.101) . □ 

Remark 7 The normal forms 3£jr + 2r]f^ and in Theorems [21 SI are valid also 

in the analytic category; see e.g. [25j. Therefore, in the analytic case, formulae (12 . 1 7ft and 
(13.9ft present Puiseux series for geodesics. 

Example 4 Consider geodesics on the surface z = y — 2x 2 in the Berwald-Moor space 
(. x,y,z ) with the metric ds = (dx dy dz)s. This yields 

F(x,y,p) — p(p — 4x), A = — 2(p 2 — Axp + 16a; 2 ), P =—4p(p + 4x), (3-17) 

and equation of geodesics (12.3ft reads 

dp = 2 p(p + Ax) 

dx p 2 — Axp + 16x 2 

The isotropic lines are solution of the differential equation p{jp — Ax) = 0. It gives two 
families of isotropic lines y = const and y = 2x 2 + const, which have the first order tangency 
on the line x — 0. Substituting them into (13.18ft . one can see that they are geodesics. 

Consider the geodesics outgoing from the point q = 0 with the double isotropic directions 
Po = 0. (Recall that for every p ^ 0 there exists a unique geodesic passing through q with 


22 



























tangential direction p, we exclude such geodesics from further consideration.) The isotropic 
geodesics y = 0 and y = 2x 2 (formula (18. 101) ) separate the (x, y) -plane into four parts: the 
upper domain y > 2x 2 , the semiplane y < 0 and two tongue-like sectors between them. See 
Fig. [4] (right), the isotropic geodesics y = 0 and y = 2x 2 are depicted with dashed lines. 

Theorem |4] claims that there exists a one-parameter family of geodesics outgoing from 
q with the double isotropic directions po = 0 into the tongue-like sectors (non-isotropic 
family (18.91) ) and there are no geodesics outgoing from q with the double isotropic directions 
p 0 = 0 into two remaining parts of the plane. Geodesics of the family (18.9|) correspond the 
admissible value v,\ — 2 (compare formulae (18.8[) . (18. 8 1) and (I3.17|) l and they can be presented 
as the Puiseux series 


y = t 6 + 3 1 6 


i>4 


di 


U-3 


P = 


dy_ 

dx 


J _dy 

3 1 2 dt 


2t 3 + dit\ where x = t 3 . 

i> 4 


Substituting the above expression for p in (I3.18|) . we obtain recurrence relations for the 
unknown coefficients a*. 

Namely, a* = 0 for all odd i (this also follows from the fact that the surface z = y — 2x 2 
is symmetric with respect to the plane x — 0). For even i we have 24a 6 + 4a| = 0, 48a§ + 
14a|a6 = 0, 72aio + 16(a4ag + a\a%) = 0, etc. In general, 

12(2i - 4)a 2 j + b 2 i = 0, i — 3,4,5,..., (3.19) 


where 62 * is a polynomial on a, 2 j with j < i with zero free term. This shows that the 
coefficient <24 is arbitrary, and all < 22 i with i > 3 are uniquely defined by equations (13.19)) . 
This gives the one-parameter non-isotropic family (13.91) . In particular, <24 = 0 gives <22 i = 0 
for all i > 3, and the corresponding solution y = x 2 presents the unique C'°°-smooth geodesic 
of non-isotropic family (13. 9 p (the corresponding value of the parameter is a = 0). 
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